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ABSTRACT 

We present VLT-ISAAC near infrared (NIR) spectro-photometric observations of 16 
post-starburst galaxies aimed at constraining the debated influence of TP-AGB stars 
on the spectral energy distribution (SED) of galaxies with stellar ages between 0.5 and 
2 Gyr, hence critical for high-rcdshift studies. Post-starburst galaxies are characterised 
by negligible on-going star formation and a SED dominated by the stellar population 
formed in a recent (< 2 Gyr) burst. By selecting post-starburst galaxies with mean 
luminosity- weighted ages between 0.5 and 1.5 Gyr and a broad range of metallicities 
(based on SDSS optical spectroscopy), we explore the parameter space over which the 
relative energy output of TP-AGB stars peaks. A key feature of the present study 
is that we target galaxies at z !=s 0.2, so that two main spectral features of TP- 
AGB stars (C-molecule band-head drops at 1.41 and 1.77 /im, blended with strong 
telluric absorption features, hence hardly observable from the ground, for targets at 
z w 0) move inside the H and K atmospheric windows and can be constrained for 
the first time to high accuracy. Our observations provide key constraints to stellar 
population synthesis models. Our main results are: i) the NIR regions around 1.41 
and 1.77 /xm (rest-frame) are featureless for all galaxies in our sample over the whole 
range of relevant ages and metallicities at variance with the Maraston (2005) "TP- 
AGB heavy" models, which exhibit marked drops there; ii) no flux boosting is observed 
in the NIR. The optical-NIR SEDs of most of our post-starburst galaxies can be 
consistently reproduced with Bruzual & Chariot (2003) models, using either simple 
stellar populations (SSP) of corresponding light-weighted ages and metallicities, or a 
more realistic burst plus an underlying old population containing up to approximately 
60% of the total stellar mass. In contrast, all combinations of this kind based on the 
Maraston (2005) models are unable to simultaneously reproduce the smoothness of 
the NIR spectra and the relatively blue optical-NIR colours in the observations. The 
data collected in this study appear to disfavour " TP-AGB heavy" models with respect 
to "TP-AGB light" ones. 

Key v^rords: galaxies: stellar content - stars: AGB and post-AGB - infrared: stars - 
infrared: galaxies - galaxies: general - galaxies: photometry 



1 INTRODUCTION 

The modelling of the spectral energy distribution (SED) of 
galaxies and its interpretation in terms of stellar popula- 
tion parameters (e.g. age, star formation history, metallic- 
ity) via stellar population synthesis (SPS) is a fundamental 
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tool to understand the properties and evolution of galaxies. 
Yet, we are far from a complete comprehension and a reli- 
able modelling of some stellar evolutionary phases which 
strongly affect the energy output of stellar populations. 
Among them, the so-called thermally pulsing - asymptotic 
giant branch (TP-AGB) phase has been the focus of debate 
among modeller s for several years (e .g. Chariot fc Bruzuall 
il991i : iMarastoiil l2005l: iBruzuall |2007| : iMarigo erall 120081 '). 
This late stage of the AGB phase for low and intermediate 
mass stars (M < 5-7 Mq) culminates in stellar populations 
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of ages between 0.5 and 1.5 Gyr. TP-AGB stars emit mainly 
in the near infrared (NIR) spectral range given the low tem- 
peratures of these stars. Quantitative predictions of the life- 
times, luminosities and spectral shapes and features of these 
stars as well as of their impact on the integrated spectra 
of even "simple" stellar populations, are still very uncer- 
tain and strongly model-dependent (e.g. iMelbourne et all 
|2012| ). One further complication derives from the fact that 
TP-AGB stars are likely to be emb edded in envelopes 
of dust produced by thems elves (e.g. iGroenewegen et all 
I2OO9I : [Ch isari fc K elson 2012), which attenuates the emerg- 
ing optical and N IR radiation by a substantial fraction 
l|Meidt et al.ll2012l ). 

The seminal work of 'Marastoiil (|2005l . Ma05 hereafter), 
who adopted the "fuel consumption theorem" approach and 
calibrated the flux contribution of this phase against optical 
and NIR photometry of Magellanic Clouds' globular clus- 
ters, dramatically pointed out how large the impact of TP- 
AGB stars on galaxy SEDs may possibly be. The Ma05 mod- 
els make two clear and tes table predictions in cont rast to 
"classical" models (such as iBruzual fc Chariot! [200^ . BC03 
hereafter), which rely on standard, less extreme prescrip- 
tions for the TP-AGB phase: in particular, i) the NIR flux 
of stellar populations of ages between «0.5 and ~ 1.5 Gyr is 
enhanced by a factor up to 3 (as a function of age and metal- 
licity); in addition, ii) strong, sharp absorption features ap- 
pear in the spectrum, especially at NIR wavelengths (and 
most notably at 1.1, 1.41 and 1.77 ^m). These features cor- 
respond to the band-heads of carbon composite molecules 
and should depend on metallicity. They represent the unde- 
batable fingerprints of TP-AGB stars as predicted by Ma05 
models. This can be clearly seen in Fig.[TJ where we confront 
the optical-NIR spectra of simple stellar populations (SSP) 
of three ages (0.5, 1 and 1.5 Gyr) at three different metallic- 
ities (0.5, 1 and 2-2.5 Zq) as derived from the Ma05 models 
(black lines) and the BC03 models (red lines) respectively. 

Observationally, detections of near-infrared CN absorp- 
tion bands have been repor ted in the integra ted spectra of 
only a few Seyfert galaxies (jRiffel et al.ll2007l ) and globular 
clusters d ominated by the stocha stic presence of bright car- 
bon stars (|Lvubenova et ai1l2012f ). In fact, from a theoretical 
point of view, carbon absorption features are not a required 
signature of AGB-star populations. Sophisticated models 
of AGB stellar evolution predict that convective dredge-up 
of carbon produced by nuclear fusion during the envelope 
thermal pulses will increase the abundance of carbon rela- 
tive to oxygen near the surface, which c an eventually lead 
to th e production of a carbon star (e.g. iMarigo fc GirardU 
[2003). While the efficiency of this process is expected to in- 
crease as the metallicity decreases, the integrated spectrum 
of most intermediate-age stellar populations could well be 
dominated oxygen-rich AGB stars without strong carbon 
absorption features. 

The largely different NIR ffux predictions for TP-AGB 
stars by different models have strong implications in terms 
of SED interpretation, leading to substantially lower stel- 
lar mass and age estimates (by up to a factor 2 or more) 
when rest-frame NIR bands are included in t he analy- 
sis o f galaxies with ages around 1 Gyr (e.g. MarastorJ 
20051: iMa raston ct al' 2006; Kanna ppan fc Gawiseij l2007l : 
Conrov et al. 2009, : Ilbc rt ct al 201(J). The popularity of 



NIR passbands as tracers of stellar mass, both at low and 



high redshift, thanks to their (generally) low sensitivity to 
age, metall icity and dust in terms of mass- to-light ratio 
(M/L) (e.g. lBeh et al.ll2003l : [zibetti et al.ll2009D . clearly calls 
for a proper assessment of the models in this respect (see 
also [Melbourne et al.|[2012l '). Verifying the reliability of SPS 
models in the regime of ages around 1 Gyr is therefore of 
paramount relevance for studies of galaxy properties and 
their evolution. 

Since 2005, a number of contrasting results and claims 
have been p roduced about the app ropriateness of the Ma05 
predictions. iMaraston et al.l (|2006t ) showed that the Ma05 
models fit the SEDs of high-z (1.5 < z < 2.5) galaxies bet- 
ter than BC03 models. However the same work shows that 
the performance of the two models differ only marginally in 
terms of goodness of fit if one allows for dust reddening, al- 
though the inferred stellar mass is substantially different. A 
better agreement of Ma05 models with the observed optical- 
NIR colours is also reported by iMacArthur et aL 

I (|2010t ). 

who used combined optical spectroscopy and optical-NIR 
photome try to constrain the sta r formation history of two 
galaxies. iLvubenova et all (|2012l ) also report detections of 
the 1.77 /im feature in two Magellanic Clouds' globular 
clusters, in agreement with Ma05 predictions. On the other 
hand, studies of Magellanic Clouds' globular clusters found 
significant discrepancies between observations and Ma05 
predictions. ILvubenova et all l|2010t ) found that the CO line 
strength (in the NIR K band) is muc h weaker than predicted 
by Ma05 in clusters of age fs 1 Gyr. IConrov fc GunnI (|2010l ) 
studied the colours of Magellanic Clouds' globular clusters 
as a function of their age and showed that Ma05 models pre- 
dict optical and NIR colours that are too red. Moreover the 
predicted age dependence of these colours does not agree 
with the data. As opposed, a substantial agreement with 
models based on much less extreme TP-AGB prescriptions 

is found. 

As already pointed out by iLancon et all (|l999l ). post- 
starburst galaxies (PSB hereafter) are the best galaxies to 
look at in order to find the fingerprints of (TP-)AGB stars 
on their SED. PSBs have a star formation history domi- 
nated by a burst in the exact age range when TP-AGB 
stars are expected to have the strongest impact, i.e. be- 
tween 0.5 and 1.5 Gyr prior to the epoch of observatiorQ. 
Their optical spectrum is characterised by a well devel- 
oped Balmer/D4000 break, the lack of strong emission lines 
(which testify a substantial cessation of star formation activ- 
ity) and by the very strong Balmer absorption lines typical of 
A stars ( hence the alternative den omination of E-l-A or K-|- A 
galaxies). IConrov fc GunnI (|2010l ) compared the optical-NIR 
colours of the sample of PSBs spectroscopically selected by 
lOuintero et al. I (I2OO4I ) to the colours of different SPS models. 
They showed (see their Fig. 16) that i — z, g — r colours are 
not reproduced by Ma05 models: not only would the blue 
portion of this colour plane be only explained by unrealisti- 
cally low metallicities, but also the age sequences in Ma05 
models pr oduce trends that are orthogonal to the observed 
age trends. iKriek et al] l|2010h performed another photomet- 
ric test, in which they fitted the stacked (rest-frame) UV- 
optical-NIR SED of 62 PSBs at 0.7 < z < 2, colour selected 



^ This holds for several different sta r formation histories (see 1 
11 and 12 in IMelbourne et al.ll2oT3) . 
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Figure 1. Comparison between SSPs at the peak of the TP-AGB phase, as derived from Ma05 models (black lines) and BC03 (red lines). 
Three different ages, 0.5, 1 and 1.5 Gyr from top to bottom, and three different metallicities, 0.5, 1 and 2 Zq (the highest metallicity 
for BC03 is actually replaced by 2.5 Zq), from left to right, are considered. Shaded regions correspond to the rest-frame spectral ranges 
covered by the H and K band ISAAC spectroscopic setups for the z ^ 0.2 galaxies presented in this work and include the strong features 
at 1.41 and 1.77 fira. 



from the NEWFIRM medium band survey (NMBS) to have 
ages around 1 Gyr, according to both BC03 and Ma05 mod- 
els. Also in this case it is found that BC03 models can re- 
produce the observed SEDs much better, in particular the 
observed blue optical- NIR colours, which are not reproduced 
by Ma05. 

In this paper we present the results of the experiment 
that we designed to test for the two main predictions of 
Ma05 models, namely the presence of sharp NIR spectral 



features and the boosted NIR flux, on a sample of PSBs 
carefully selected from optical spectroscopy and newly ob- 
served in NIR spectroscopy (H and K band) with ISAAC at 
the ESO-VLT. The selection criteria and properties of the 
sample are fully detailed in Section [S] Here we would like 
to stress the uniqueness of this work, i) For the first time 
we cover both the optical and the NIR range of PSBs with 
flux calibrated spectra to look directly for the NIR spectral 
features at 1.41 and 1.77 ^m predicted by the Ma05 mod- 



© 2012 RAS, MNRAS OOO.mfSTI 



4 S. Zibetti et al. 



els and resulting from the large impact of TP-AGB C-rich 
stars. These features are unaccessible or extremely difficult 
to measure from the ground as they roughly coincide with 
the atmospheric gaps between J and H, and between H and 
K bands, respectively, for nearby objects. Therefore we ob- 
served PSBs in the redshift range 0.15-0.25 so that the C- 
features move well into the observable windows and are min- 
imally affected by uncertain telluric corrections, as shown in 
Fig. HI a) Contrary to previ ous works based either on colour 
selection (jKriek et al.ll201C ) or on a purely phen omenologi- 
cal spectroscopic selection ( Conrov fc Gunnl20ld l. the selec- 
tion of our sample is based on both phenomenological crite- 
ria (h igh Balmer line strength and low emission lines, iGotd 
[2003) and bayesian estimat es of light weighted ages from 
optical spectral indices as in iGallaz zi et a ir (!2005f ). This al- 
lows us to focus our observations and analysis on the precise 
range of stellar ages where TP-AGB stars are expected to 
contribute the most, for a variety of stellar metallicities. 

The paper is organised as follows. In Section [2] we 
present the sample, the new ISAAC NIR spectroscopic ob- 
servations and the data reduction, including spectrophoto- 
metric calibrations and integration with optical data. In Sec- 
tion [3] we present the results and discuss possible spurious 
effects and contaminations that might affect our conclusions. 
A comparison with the analysis of Krick ct al. (2010.) is also 
presented. Finally, in Section 3] we summarise our findings 
and propose our conclusions and future developments. 



2 SAMPLE AND OBSERVATIONS 

2.1 Sample selection and characterisation 

We draw our PSB targets from the spe ctroscopic sample 
of the Sloan Digit al Sky Survey (SDSS, lYork et al.l I2OO0I: 
Strauss et al.l |2002| ) data release 7 (DR7. [Abazaiian et al.l 



20091 ) ■ The initial pre-selection was b ased on the l ine equiv- 



alent width (EW) criteria defined in ICotol (|2005h . namely: 
EW(m)> 5 A, EW(Ha)> -3 A, EW([OII])> -2.5 A0 In 
order to ensure that the two strong NIR features correspond- 
ing to the band-heads of CN (1. 41 /jm) and C2 (1-7J Atm) 
due to carbon rich AGB stars (cf. lLancon fc Mouhcinell2002l ) 
are observable from the ground, we require the galaxies to be 
at 0.15 < z < 0.25, so that 1.41 /xm and 1.77 fira (at rest) 
move into the H- and K-band windows, respectively. This 
selection is based on the SDSS-DR7 casjobs catalog We 
further restrict the sample by applying the visibility cuts 
in order for targets to be observable from Paranal during 
the winter semester: DEC< 20 deg and RA between and 
230 deg or RA> 330 deg. Finally we require that the stellar 
light-weighted age of the galaxies is in the range 0.5 to 1.5 
Gyr, i.e. at the maximum of the predicted TP-AGB lumi- 
nosity contribution, based on the estimates computed as in 
IGallazzi et al.l (|2005l ') and available in the MPA-JHU cata- 
lo£|. Given the relative rarity of the post starburst phase and 



^ We adopt the convention of positive EW for absorption and 
negative EW for emission. For the [O II] emission, the sum of the 
doublet is considered. 
^ http://casjobs.sdss.org7' 

* http://www.mpa-garching.mpg.do/SDSS/ Note that although 
the MPA-JHU catalog is limited to DR4 only, no loss of possible 
targets derived from this cut. 



the restrictive age, coordinate and redshift cuts (the latter 
particularly relevant as it is almost at the extreme bound- 
ary of the SDSS main galaxy sample), the final sample is 
composed of sixteen galaxies only. 

The properties of the sample are summarised in Table 
[T]and in Fig. [21 This figure illustrates the distribution of our 
galaxies in the classical HcSa vs D4000n plane. In this plane 
the main sequence of galaxies with smooth continuum star 
formation history is displayed by the grey intensity levels, 
which represent the distribution of SDSS-DR4 galaxies of 
increasing ages from top-left to bottom right. Galaxies ly- 
ing above this sequence, such as our galaxies, represented 
by the red dots with error bars, are characterised by rela- 
tively recent bursts of star formation with ages between a 
few hundred Myr to a couple of Gyr, which dominate the 
light of the galaxy. For comparison, tracks of Simple Stellar 
Populations (SSP) with different metallicities from BC03 are 
shown; the typical bell-shape of these SSP tracks peaks in 
II(Sa at ages of 0.3-0.5 Gyr (depending on metallicity) and 
decreases to H^a ~ 3-5 for ages of a few Gyr. By design, 
our sample populates the intermediate range of burst ages 
in this diagram. 

The statistical estimates of light-weighted stellar ages, 
metallicity and mass computed from optical spectral absorp- 
tion indices (D4000„, RSa - j-HiA, H/3, [MgFe]', [MgaFe]) 
with the bayesian method of lGallazzi et al.l ( 20051 ) and used 
in the selection are reported in Tab. [T] (columns 6, 7 and 8), 
along with the ID (1), coordinates (2, 3), redshift (4) and 
apparent r-band petrosian magnitudes (5) of the galaxies 
from the SDSS. The sample comprises intermediate mass 
galaxies, between ~ 3 and ~ 7 x lO^'^Mo, with ages in the 
range 0.8 — 1.4 Gyr and metallicities between 0.2 and 2.3 so- 
lar. Most galaxies have indeed Z > 0.5Zq, as expected given 
their stellar mass. Despite the paucity of galaxies at low age 
(< 1 Gyr) and low metallicity (< IZq), the sample provides 
a good coverage of the relevant ages for the TP-AGB phase 
and of the metallicities most typical for galaxies around the 
characteristic luminosity L*. 

We note that the estimates of the stellar population 
parameters of IGallazzi et al.l (|2005l ) rely on the compari- 
son with synthetic spectral libraries based on BC03 mod- 
els, therefore they could be biased with respect to esti- 
mates based on other synthesis models, most notably the 
Ma05 ones. However, optical spectral indices short-wards 
of 5500 A are generally well agreed upon among different 
modellers and, most importantly, are only weakly affected 
by TP-AGB stars, which mainly influence the NIR spectral 
region. Thus age and metallicity estimates from optical ab- 
sorption features are robust against different TP-AGB pre- 
scriptions and the use of different models. This is even more 
the case for spectra dominated by A stars, where Balmer 
absorptions are very strong. Furthermore, in Sec. |3| (|3.2l and 
I3.3l in particular) we thoroughly explore if and how different 
(combinations of) ages and metallicities can possibly explain 
the observed properties of our galaxies, irrespective of their 
original selection. 

Based on the SDSS optical images and, to minor ex- 
tent, on our ISAAC NIR images, we note that the ma- 
jority of galaxies display smooth and centrally concen- 
trated light distributions, typical of early type galaxies, with 
the only exception of PSB J0227— 0015, which has an ex- 
tended and irregular "halo". However, other four galax- 
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Table 1. Sample properties. 



ID RA Dec z r'potro L-w Age L-w Z 
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Figure 2. Sample properties in the classical H(5-D4000n diagram, 
a proxy for burstiness vs. age of the stellar population. Red points 
with error bars are the sixteen galaxies selected for the present 
study; the underlying grey-scale image represents the density dis- 
tribution of a complete sample of galaxies from SDSS-DR4. As 
one can see, our sample selects galaxies dominated by a burst in 
a range of young to intermediate ages. Note that the estimates of 
H5a used in this plot are not exactly the same as the E.W. mea- 
surements used in the sample selection (therefore values < 5 A are 
not inconsistent with the selection criteria reported in the text). 



ies (PSBJ1039+0604, PSB J1046+0714, PSBJ1119+1313 
and PSB J1206+0918) have close and possibly interacting 
companions; in many cases, faint structures (like shells or 
streamers) can be tentatively seen. This suggests a relatively 
recent merger or interaction as a possible origin of the ceased 
starburst in these galaxies. 

2.2 Observations and data reduction 

The sixteen galaxies of the sample were observed in the pro- 
gramme ID 086.B-0733 (P.I.: S. Zibetti) at the UTS Melipal 
of the ESO- VLT at the Paranal O bservatory. The ISAAC 
instrument l|Moorwood et all 1 19981 ) has been used in two 
different spectroscopic setups to obtain low resolution spec- 
tra in the H and K bands. 1.5 arcsec slit and SWSl-LR grism 
were used in both setups, with central wavelengths of 1.65 
and 2.2 /im respectively. This gives a spectral resolution of 
f» 350 (50 A) and fa 300 (70 A) in H and K respectively. Ob- 
servations were carried out in service mode with clear (not 
necessarily photometric) sky and maximum seeing of 1.4 arc- 
sec to minimise slit losses and match the fibre size and see- 
ing of the SDSS optical spectra as well as possible. Standard 
ABBA on-slit nodding with 70 arcsec offsets and 20 arcsec 
jittering were adopted to allow accurate sky subtraction and 
bad pixels and cosmic ray rejection. Eight exposures of 180 
sec each were taken in H band and six exposures of 150 sec 
each were taken in K, for total 24 minutes of integration in H 
and 15 minutes in K per target. This eventually yields typi- 
cal signal-to-noise ratio (SNR) of « 50 and 25 per resolution 
element in H and K respectively. Spectro-photometric stan- 
dard stars were observed with identical setup immediately 
after or before the science observations to provide a reliable 
estimate of the spectral sensitivity function and compute the 
corrections for telluric absorption. Standard data reduction 
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procedures have been implemented in IRAE0. Wavelengtii 
calibration has been performed based on the sky OH lines 
as measured on a median combined image of the scientific 
frames. This approach turned out to be much more stable 
and reliable than using the arc lamp calibrations taken in 
day time. We have tried different approaches to optimally 
correct the spectra for telluric absorption. As our best solu- 
tion we adopted the following strategy: the empirical high- 
resolution transmissivity function of the average atmosphere 
above Paranal (provided by ESO) is finely tweaked by scal- 
ing and changing the resolution until we can exactly com- 
pensate for the observed absorption at various wavelengths 
on the spectrum of the spectrophotometric standard; the 
coadded spectrum of each galaxy is then corrected for the 
same absorption pattern as for the closest (in time and air- 
mass) standard star. We have applied the telluric correction 
to both standard stars and galaxies before computing the 
sensitivity function and applying the flux calibration. 

In order to assemble a SED including the visible range 
and the H and K NIR bands, absolute flux calibration is re- 
quired. To this goal we have obtained high SNR images in H 
and Ks bands in photometric conditions with ISAAC as part 
of the same program. "Total" integrated magnitude derived 
from these images are used to rescale the spectra. By "to- 
tal" magnitude we adopt the magnitude integrated within 
the petrosian aperture, as defined by the SDSS photometric 
pipeline on the SP SS r-band images in data release 8 (DR8, 
lAihara et al.llioill 'FI. DR8 is chosen as it fixes a number of 
issues especially with sky subtraction. H and Kg magnitudes 
are used to rescale the respective ISAAC spectra in the NIR, 
while the optical spectrum from the SDSS is normalised ac- 
cording to the r-band petrosian magnitud^. Conversions 
between NIR magnitudes natively calibrated in the Vega 
system for the ISAAC filters into the corresponding flux 
densitie^ for o ur spectra are computed based on the spec- 
trum of Vega of lKurucj (|l992l '). distributed by BC03. Typical 
total photometric uncertainties for ISAAC petrosian magni- 
tudes are estimated ~ 0.05 mag and include photometric 
zero point fluctuations (~ 0.03 mag, derived from repeated 
observations of standard stars during the nights) and back- 
ground subtraction uncertainties For the SDSS we assume 
errors of 0.02 mag, which is the maximum reported for our 



^ IRAF is distributed by tlie National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities 
for Researcii in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 

^ Data are obtained from the DR8 casjobs catalog server 
[http: / /skyser vice .pha.jhu.edu /casj obs 

' We adopt the petrosian magnitudes distributed in the SDSS- 
DR8. We have independently checked these values by doing 
our own aperture photometry on the SDSS images and found 
substantial agreement: our estimates are 0.011 mag fainter on 
average, with a scatter of 0.015 mag. Only for two galaxies 
(FSB J1133-I-0205 and PSB J0151-0056) we find fainter magni- 
tude by more than 0.03 mag (-1-0.040 and -1-0.033 respectively). 
These small discrepancies are well within the errors adopted in 
Sec. 132] and [33] 

* Average fiux densities within broad band filters are com- 
pu ted according to the integral defining AB magnitudes given 
bv lFukugita et al.1 l| l99d ). 

^ Photom etric measurements similarly performed on the 
UKIDSS jLawrence et al.]|2007l ) images that are available for a 



objects by SDSS-DR8. Foreground Galactic extinction cor- 
rection is finally appl ied using the optical depth derived from 
ISchlegel et ahl (1 19981) and assumi ng the standard e xtinction 
laws of ICardelh et al.l (| 19891 ) and lO'DonneUl (|l994l ). 



2.3 Aperture effects 

The spectra of our galaxies, both SDSS and ISAAC, are 
taken in relatively small apertures or slits centred on the 
brightness peak: 3-arcsec fibres are used in the SDSS with a 
typical seeing of ~ 2 arcsec while 1.5-arcsec long-slits with 
seeing of 1.4 arcsec or better are used with ISAAC. Rela- 
tive aperture size and seeing are designed to conspire and 
collect a similar fraction of the galaxy light. By rescaling 
the spectra according to the "total" petrosian magnitudes 
we implicitly assume that stellar population gradients are of 
limited relevance and the ce ntral region o f the galaxies are 
representative of the whole. jPracv et al.l (|2012l ) have anal- 
ysed a small sample of post-starburst galaxies with integral 
field spectroscopy and showed that, in reality, a large frac- 
tion of them have cores that are younger than the rest of the 
galaxy, hence our assumption does not strictly hold. Unfor- 
tunately, because of different seeing conditions and geome- 
try of the apertures, a precise and homogeneous rescaling of 
the spectra is very hard to achieve. However, by comparing 
fiuxes and colours that are obtained in the petrosian and in 
the fibre apertures we can show that the aperture effects for 
our sample are small enough not to affect our conclusions. 

First of all, we note that at the relatively high redshift 
of our sample the apparent size of the galaxies is typically 
2.5 arcsec (median petrosian radius, maximum 3.1 arcsec), 
so that the 3-arcsec SDSS fibre typically collects already 
half of the total flux. The ISAAC spectra roughly collect 
the same fraction of flux. The long slit geometry in this case 
integrates into the spectra not only the central part, but also 
the outskirts of the galaxy, so that the ISAAC spectra are 
indeed even more representative of the whole galaxy. 

In second place, we have checked by how much "total" 
petrosian colours differ from flbre colours in the optical. Fi- 
bre colours are computed in such a way that they correspond 
as closely as possible to the SED portion effectively captured 
in the SDSS spectra. The strongest effects are expected and 
observed in g — i, since this colour roughly corresponds to 
u — r aX z ~ 0.2 and therefore is the most sensitive to stel- 
lar population properties by bracketing the D4000 break. 
We flnd that, on average, fibre g — i is bluer by 0.027 mag 
with respect to the corresponding total petrosian colours, 
with a rms of 0.037. Only four galaxies (out of 16) have 
a fibre colour which is bluer than total by more than 0.05 
mag: PSB J1046-h0714, PSB J1015+0103, PSB J0227-0015 
and PSB J0328-I-0045. We dub these galaxies "blue core" 
in the following. For the remainder, differences in colour 
can be considered as not significant, although the data hint 
at starbursts preferentially occurring in the central regions 
of galaxies r ather than in the outskirts, consistently with 
iPracv et al.l (|2012l ). Looking at r — i, which roughly corre- 
sponds to the rest-frame optical colour mo. 55 — mo. 70 that 



subsample of our galaxies are found in good agreement within 
the quoted uncertainties. 
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will be discussed in the following, we find even smaller dif- 
ferences of 0.012 mag on average (0.021 mag rms). These 
are well within the spectrophotometric calibration errors. 

Finally, we estimate the relative aperture corrections 
between the optical and NIR bands. Specifically, we con- 
sider the aperture corrections in r and //-band, defined by 
ap_corr = mflber — "T-petro- For the r-band we use the quan- 
tities provided by the SDSS database. For the //-band we 
use the total petrosian magnitude as derived above, while 
to obtain fibre magnitudes we first convolve the ISAAC 
images to a common PSF of FWHM=2 arcsec and then 
integrate the flux within a circular aperture of 3 arcsec- 
ond in diameter, as done in the SDSS pipeline. Since for 
most of the ISAAC images we lack a sufficient number 
of high SNR stars to compute an accurate PSF, the final 
FWHM is approximated within 0.1 arcsec. We find an av- 
erage difference of aperture corrections in r and //-band 
< ap_corrr — ap^corrn >= 0.026 mag with rms of 0.1 mag, 
indicating that on average our galaxies are slightly more 
concentrated in r than in H, although with very marginal 
significance. This is consistent with the existence of few blue 
cores and, in general, of very small systematic SED varia- 
tions between the spectral and the total apertures. 

This analysis of the aperture effects justifies our choice 
to use the total petrosian magnitude to compute the rela- 
tive normalisation of the optical and NIR spectra. Total pet- 
rosian magnitudes are much more stable than magnitudes 
integrated in small apertures comparable with the seeing, 
which therefore rely on uncertain PSF convolution. On the 
other hand, we have shown that they do not introduce sig- 
nificant biases between different spectral regions. The sys- 
tematic aperture bias of 0.026 mag between r and H is in 
fact well within the photometric uncertainties estimated in 
the previous section. 

In the following analysis we will adopt 0.1 mag as stan- 
dard error on the optical-NIR colours and 0.05 mag on the 
optical colours, to take into account both photometric errors 
and aperture effects. 



3 RESULTS AND ANALYSIS 

3.1 Composite optical-NIR spectra: absence of 
sharp NIR features 

The composite optical-NIR spectra obtained as explained 
in the previous section are presented in Fig. [3l in the rest 
frame and normalised to the flux density at 5500 A, in or- 
der to ease a prompt comparison between all galaxies. The 
actual normalisation refers to the median flux in the range 
5480-5520 A, in order to minimise the effect of noise on the 
individual spectral pixels. The SDSS portion of the spec- 
tra covers approximately the rest-frame range 3000-7500 A, 
ISAAC H 1.25-1.50 fim, and ISAAC K 1.65-2.00 fim. For 
each galaxy, three panels display the full optical-NIR spec- 
tral range (top left plot), a zoom into the optical region 
that includes the most significant diagnostic lines (top right 
plot), and the NIR H-K range (bottom plot). Green lines 
are used to plot the spectra at high resolution: the original 
Ri 40 and 50 A resolution at ~ 4 and 6 A per pixel in 
the rest frame is used for the H and K spectra respectively, 
while for the SDSS spectra (original resolution 2.5 A at 



« 1 A per pixel in the rest frame) a smoothing over 3 pix- 
els is applied to make the plot more readable. Black thick 
lines are used to (over)plot heavily smoothed versions of the 
spectra, in order to ease the identification of possible spec- 
tral features, especially in the NIR bands: smoothing over 
15 pixels is adopted both in the optical and the NIR spectra, 
resulting in rest frame resolution of f» 15, 60 and 90 A in 
optical, H and K ranges, respectively. Vertical dotted lines 
mark the position of the main (expected) spectral features: 
the Balmer series (up to Hr; 9-2) and the two features of CN 

1 1.41 /xm) and C2 (1-77 /^m ) due to carbon rich AGB stars 
Lancon fc Mouhcinell200^ ). For each galaxy in the bottom 
panel we report also the light- weighted stellar age and metal- 
licity and corresponding u ncertainties derived from optical 
absorption features as in Ga llazzi et al., (20051 . 

All galaxies display typical E-l-A or post-starburst spec- 
tra with very strong Balmer absorptions, well developed 
D4000 break and blue continuum. The actual slope of 
the optical continuum generally shows the well known de- 
pendence on age and metallicity, although dust attenu- 
ation may affect the slope as well. We do not attempt 
to put any constraint on this parameter, but discuss its 
possible impact on our conclusions below. Three galaxies, 
PSB J0227-0015, PSB J0328-0045, and PSB J1150-0107 
show significant Ha emission (and some higher ionisation 
lines for PSB J0328— 0045), indicating some residual ongo- 
ing star formation, which is also confirmed by (marginal) 
detections at 22 fim (W4) in the WISE AU Sky Survey 
(|Wright et al.l l2010'l. measured with SNR of 4.4 and 2.8, 
and 2.0, respecti vely. Howeve r, the equivalent width in all 
cases matches the lCotol l|2005h criteria, being — 2.27±0.10 A, 
-2.80 ± 0.10 A, and -2.43 ± 0.24 A, respectively, aU larger 
than —3 A. Given the observed strength of all Balmer ab- 
sorption lines, the contribution of newly formed stars (or 
AGN) can be largely neglected in comparison to the radia- 
tion by the intermediate-age population we aim at study- 
ing here. It is worth noting that two of these galaxies 
(PSB J0227-0015 and PSB J0328-0045) also have a blue 
core, supporting the idea that blue cores might be related to 
some residual star formation extending close to the present 
time. The absence of emission lines and blue cores from the 
majority of the galaxies in our sample further confirms their 
nature of genuine post-starburst galaxies. 

What is most important for the scope of this work 
is the lack of any evidence for sharp features in the 
NIR in any of the 16 galaxies. Most of them display 
completely featureless continua with approximately uni- 
form slope within the two NIR windows. Some tentative 
evidence for a wavelength-dependent slope is seen in a 
few objects: PSB J1039-h0604, PSB J1046-h0714 and espe- 
cially PSB J1314-I-0430 appear to change their slope around 
1.41 fim, while in PSB J1119-H313 and PSB J1133-H0205 a 
change of slope can be tentatively seen around 1.77 /im. 
Similar mild changes of slope, most likely due to molecular 
bands in the atmospheres of cool stars, are compatible with 
BC03 but are in fact substantially dissimilar from the sharp 
drops by a few 10% in flux predicted by Ma05 models (see 

Fig. DO). 

To further support these conclusions, we stacked all 16 
spectra in H and K separately, after normalising them in a 
1000 A-wide region around 1.41 and 1.77 fim, respectively. 
The results are displayed in Fig.|4l with the black lines show- 
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Figure 3. Composite rest-frame spectra for the sixteen galaxies of the sample, including optical SDSS spectra and the new NIR H and 
K band spectra observed with ISAAC at the ESO-VLT. Correction for foreground Galactic extinction and normalisation at 5500 A are 
applied after relative normalisation between the optical and the NIR spectra based on SDSS r-band photometry and our own ISAAC 
imaging, as detailed in the text. The full/high resolution version of the spectra is plotted in green (rest-frame resolution 2.5 A for the 
SDSS and 40 and 50 A for ISAAC); the thick black curves are smoothed versions of the same spectra, to filter out the noise and improve 
the visibility of the features. The resulting effective resolution in these smoothed spectra is 15, 60 and 90 A rest-frame in optical, H 
and K ranges, respectively. The new NIR H and K band spectra observed with ISAAC at the ESO-VLT are displayed in the bottom plot 
of each panel. The full optical-NIR range is shown in the top-left plot, while a zoom-in of the optical range is shown in the top-right plot. 
Vertical dotted lines mark the position of the main (expected) spec tral features: the Balmer s eries (up to Hr; 9-2) and the two features 
of CN (1.41 fj,m) and C2 (1.77 fim) due to carbon rich AGB stars llLancon &: Mouhcine|[2"002l 'l. Light-weighted ages and metallicity are 
reported for each galaxy (sec text for details). 
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Figure 3. (continued) 



ing the average spectrum and the shaded green area indicat- 
ing the 16-84% inter-percentile range (approximately ±lcr 
range). This figure shows that on average carbon molecules 
in the atmospheres of evolved stars produce band-head drops 
that are definitely less than 5% of the continuum. 

3.2 Optical-NIR colours vs optical absorption 
indices 

Stellar absorption indices in the optical range are ex- 
tremely powerful diagnostics of stellar population proper- 
ties. Decades of development and tests have led to sub- 
stantial agreement on their interpretation in terms of ages 



and global metallicity. In this section we exploit part of the 
wealth of information that we can extract from the opti- 
cal SDSS spectra of our galaxies for two main purposes: i) 
probe the ability of "TP-AGB light" models (represented 
by BC03) and "TP-AGB heavy" models (Ma05) to repro- 
duce the observed trends and distributions in the parameter 
space of optical absorption indices and optical-NIR colours; 
ii) investigate whether the lack of NIR spectral features can 
be due to selection biases such that we misinterpreted op- 
tical diagnostics (e.g. due to the use of BC03-based model 
libraries) and ended up missing the relevant phase for TP- 
AGB stars. 

We start considering only age-sensitive absorption in- 
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Figure 3. (continued) 



dices and optical-NIR colours and build the diagnostic plots 
presented in Fig. [5] for BC03 and Ma05 in the top row and 
bottom row, respectively: coloured symbols and lines show 
the predicted optical-NIR colour mo. 55 — mi. 40 as a function 
of H(Sa + H7A, D4000n, and the optical colour mo. 55 — mo. 70 
(from left to the right ) r'^i For BC03, ages range from 50 Myr 
(lower ages are not relevant for the present analysis as they 
would yield much bluer colours than observed) to 13.6 Gyr. 
Metallicities range from 1/50 solar to 2.5 times solar, and are 



Here mo. 55, rno.70 and mi. 40 denote the AB magnitudes corre- 
sponding to the flux density around 0.55, 0.70 and 1.40 iim (rest), 
respectively. 



represented with different symbols and colours from blue to 
red (solar in orange, see figure caption). Symbols of increas- 
ing size are used to mark discrete increasing ages, with the 
crucial age of 1 Gyr identified by a star. The spectral indices 
H(5a -I- H7A and D4000n are measured after convolving the 
model spectra to match the resolution and average velocity 
dispersion of the galaxies (< a >= 164km s~^). The Ma05 
SSP spectra are distributed over the full optical-NIR range 
only at low resolution (20 A sampling), therefore cannot 
be use d to compute reliable in d ices, H^a -I- H7A in partic- 
ular iMasiiioZ^IliliffiEicS (|201ll ) delivered high reso- 

The effect of such a coarse sampling is to underestimate 
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Figure 3. (continued) 



lution spectra for a subsample of the Ma05 age-metallicity 
grid, over the optical range: only three metallicities are avail- 
able (0.5, 1.0 and 2.0 solar), with the lowest ages being 200, 
30, 400 Myr, respectively. We use these models (in their 
STELIB version), convolved to the effective resolution of 
our data, to compute the H(5a + H7A and D4000n plotted in 
the first two panels of the bottom row of Fig. O optical-NIR 
colours are computed from the original Ma05 SEDs of cor- 
responding age and metallicity. For the colour-colour plot 



H5a + H7A by up to 
to ^ 0.1. 



5 A and to overestimate D4000n by up 



the full grid of Ma05 is used instead, which allows extend- 
ing the metallicity down to 1/20 solar and the age down 
to 50 Myr as for BC03. Observational data points derived 
from our composite optical-NIR spectra are plotted as error 
barsH 



Errors on the colours are uniformly set equal to 0.05 and 0.1 
mag for the optical colour mo. 55 — mo. 70 and the optical-NIR 
colour mo. 55 — mi. 40 respectively, following the discussion in Sec. 
12.31 Errors on indices are derived from the MPA-JHU catalog and, 
in the case of D4000n, incremented a posteriori to take into ac- 
count flux calibration errors, as statistically derived by comparing 
duplicate observations of a large sample of SDSS galaxies. 
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Figure 4. Stacked H {left panel) and K (right panel) spectra of the 16 galaxies. Each spectrum has been re-normaUzed to 1 in a 
wavelength range of 1000 A around the expected positions of the features at 1.41 and 1.77 /xm (marked with red dashed vertical lines), 
respectively. Band-head drops of more than a few per cents on average can be ruled with very high confidence. 



As shown by the top row of Fig. [S] BC03 SSPs provide 
a good match to the observed spectral diagnostics, with the 
only exception of a couple of galaxies, PSB J1150— 0107 and 
PSB J1206-I-0918, which appear to have a too blue mo. 55 — 
mi. 40 NIR-optical colour for the given mo. 55 — mo. 70 optical 
colour (or, reversely, a too red mo. 55 — mo. 70 for the given 
mo. 55 —mi. 4o) and will be discussed in Sec. 13.41 In particular, 
galaxies cover the range in age and metallicity corresponding 
to light-weighted ages and metallicities listed in Table [U 
which are mostly relevant for the TP-AGB phase. 

As far as the comparison with the Ma05 models is con- 
cerned, we note that the observed age-sensitive indices ap- 
proximately match the range of ages given by our light- 
weighted estimates. Most of the observed H5a + H7A and 
D4000n values correspond to SSP ages between 0.5 and 1 
(1.5) Gyr for the solar (sub-solar) metallicity track, while 
ages smaller than 0.5 Gyr would be inferred from the 2 
times solar metallicity track. For the given indices, how- 
ever, the predicted mo. 55 — mi. 40 NIR-optical colours are 
on average 0.5 mag redder than the observed ones. This is 
consistent with the Ma05 "TP-AGB heavy" models having 
the NIR very much boosted with respect to BC03 in the 
age range around 1 Gyr implied by the Balmer indices. In 
the rightmost panel of Fig. [S] bottom row, we see that, un- 
less a metallicity as low as 1/20 solar is representative of 
most galaxies and ages are significantly larger than 1 Gyr, 
the Ma05 models predict mo. 55 — mi. 40 NIR-optical colours 
that are on average 0.5 mag redder than observed, for the 
observed optical mo. 55 — mo. 70 colours 

We note that the inclusion of dust with a standard ef- 
fective attenuation curve in the models does not alleviate 
the tension between the Ma05 models and observations. We 



Similar conclusion is reported bv lConrov fc GunrJ l|201Clh con- 
cerning the colours of their sample of PSB galaxies. 



illustrate this in Fig. [5] by drawing arrows which represent 
the offset to apply to each model when an optical extinction 
of Av ~ 1 ma g and an effectiv e attenuation law t\ oc A~°'^ 
(foUowing .Chariot fc Fallll2000l . for an SSP) are applied. As 
expected, the effect of dust is to move model tracks to redder 
colours: this would worsen the disagreement in the index- 
colour planes. The reddening vector is approximately paral- 
lel to the model tracks in the colour-colour plane: accounting 
for dust attenuation in the models would not substantially 
change the (dis)agreement with the observations in terms of 
colours, although it would require even lower metallicities 
and ages. 

The hypothesis of a metallicity as low as 1/20 solar is 
not only unlikely to apply to a full sample of relatively mas- 
sive galaxies with M* « 5 x lO^^M©, but is also inconsistent 
with direct metallicity constraints which can be derived from 
absorption indices. In Table[T]we already showed metallicity 
estimates derived fro m the vast library (ba sed on BC03) and 
bayesian method of iGallazzi et al] l|2005l ): these estimates 
range from 0.2 Zq minimum to 2.3 Zq maximum, roughly 
uniformly distributed. In Fig. |6]we show that metallicities as 
low as 1 /20 solar are indeed ruled out already from a simple 
comparison of standard absorption line indices on which dif- 
ferent models are in substantial agreement, namely [MgFe]' 
and H/3. [MgFe]' in part icular is chose n as o ptimal metal- 
licity indicator following iThomas et al.l l|2003l ) who demon- 
strated its insensitivity to a/Fe abundance ratio; H/3 serves 
to constrain the age of the stellar population and hence 
alleviate the well known age- metallicity degeneracy. Fig. [6] 
shows the indices measured in our galaxies as crosses with 
error bars, while tracks of different colours are derived from 
BC03 SSPs in the age range 50 Myr - 13.6 Gyr (symbols 
are the same as in Fig. [5}. Four different metallicities, 1/50, 
0.4, 1 and 2.5 solar are plotted in purple, blue, green or- 
ange and red, respectively. It is apparent from this plot that 
metallicities significantly lower than 0.4 solar are inconsis- 
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Figure 5. Comparison between tlie observed properties of the post-starburst galaxies in our sample (points with error bars) and a suite 
of SSP models (coloured symbols and lines) from BC03 (top row) and Ma05 {bottom row). The optical-NIR colour mo. 55 — mi. 40 is 
plotted against H<5a.+H7a {left panels), D4000n {central panels), and the optical colour mo. 55 — mo.70 {right panels). Each track (blue 
with upside-down triangles, green with diamonds, orange with circles and red with crosses) corresponds to metallicity of 1/50, 0.4, 1 and 
2.5 solar, respectively, for the BC03 models in top row, and to metallicity of 1/20, 0.5, 1 and 2 solar, respectively, for the Ma05 models 
in bottom row. The tracks for BC03 and the colour-colour plot of Ma05 extend from 50 Myr to 13.6 Gyr. As Ma05 only provides low- 
resolution SEDs, H5a-I-H7a and D4000n for these models are taken from the corresponding STELIB-based SSPs published by Maraston 
& Stromback (2011), limited to three metallicities, 0.5, 1.0 and 2.0 solar, and starting with ages of 200, 30, 400 Myr, respectively. Stars 
mark 1 Gyr, while symbols of increasing size mark ages of 0.3, 0.5, 1.5, 2 and 10 Gyr. IMF is assumed Salpeter in all cases. The arrow in 
each plot shows the shift produced by Ay = 1 mag for a uniform attenuation following tx oc A""'^. While colours and spectral indices 
can be reproduced simultaneously by BC03 SSPs, Ma05 SSPs yield optical-NIR colours that are inconsistent with age-sensitive spectral 
indices and optical colours, unless metallicities as low as 1/20 solar are advocated for all galaxies. 



tent with the data-points. We have checked that the same 
conclusions are obtained if indices based on BC03 models 
are rep laced with those computed by Thomas an d collabo- 
rators (jThomas et al.ll2003l . [20o3 : lKom et"al]l2005l ): the only 
difference between the two sets of models is that Thomas' 
models shift the lower limit allowed by the data to slightly 
higher metallicities. 

3.3 Effects of composite stellar populations 

A possible explanation why the strong TP-AGB features 
predicted by Ma05 models are not seen in our spectra is 
that in the NIR there is a substantial contribution by old 
stars, which are almost completely outshined by the younger 



population in the optical bands. In fact, it is quite unlikely 
that the entire stellar population of galaxies as massive as 
several 10^°Mq is formed in a burst ~ 1 Gyr old and a 
more realistic scenario is that the burst shines atop an old 
component. To explore this hypothesis we have repeated the 
test presented in Fig. [5] using composite stellar populations 
(CSP) instead of SSPs. 

In a first experiment, each CSP is composed of two com- 
ponents: a fixed 10-Gyr old, solar metallicity SSP (which 
we will call "old component" in the following) and an SSP 
of variable age, from 50 Myr to 13.6 Gyr, and metallicity, 
from 1/20 to 2-2.5 solar (which we will call "young" com- 
ponent or "burst", and corresponds to the SSPs of Figure 
[Sjl. In the following we refer to these models as CSP-S, 
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Figure 7. Same as FigureO but now using tlie composite SPs tliat are referred to as CSP-S in the text: a 10 Gyr old, solar metallicity 
component is combined with SSPs with the same properties as (and represented as) in Figure [5] The ratio of masses as formed in the 
fixed "old" vs variable "young" components is 9:1. For both BC03 and Ma05 models, such a combination can reproduce the spectral 
indices with "young" components younger than ^ 0.5 — 1 Gyr (even more so for the Ma05 models); however the optical-NIR colours 
cannot be reproduced by these models. Inclusion of standard dust attenuation (illustrated by the arrows) does not help also in this case. 



where "S" indicates that the metallicity of the old com- 
ponent is solar. This choice appears reasonable consider- 
ing that Z > Zq is typical for the old component of mas- 
sive Ralaxies with M * > 1O^°M0 (e.g. iGallazzi et al-llioosl : 
omas et al.ll2005l ). We have varied the relative fraction of 
(initial) mass in the old component and in the burst from 0:1 
to 0.9:0.1. Fig. [7|shows the same diagnostic plots we already 
used to analyse pure SSPs (Fig. [5]), but now plotting CSP-S 
models for the extreme "contamination" , i.e. for 90% of the 
formed mass in the ol d component and 10% in the burst (as 
was also assumed by IConrov fc GunnI |2010| ) . Symbols and 
colours are the same as in Fig. [S] and refer to age and metal- 
licity of the variable burst component alone. CSP-S models, 
either BC03 or Ma05, can easily reproduce the distribution 
of data-points in H5a -|- H7A and D4000n: the contamina- 
tion by the old component is compensated by a younger 
burst component (< 500 Myr) with respect to those that 
provided the best match in the pure burst (SSP) models. 
Such a combination would easily explain why we do not see 
the sharp NIR features due to TP-AGB stars predicted by 



Ma05: not only part of the NIR emission is provided by the 
old and TP-AGB free stellar population, but also the young 
component is younger that the peak of the TP-AGB phase. 

However, problems arise when optical-NIR colours are 
considered. The addition of an old component at solar metal- 
licity to a young burst yields significantly redder optical- 
NIR colours with respect to a pure burst. As a consequence, 
BC03 models, which provided a decent match in the case 
of a pure SSP, are now shifted to inconsistently red colours 
with respect to the observations. In the case of Ma05 mod- 
els, which already had troubles being too red in optical-NIR 
colours at almost any age, they become inconsistent with the 
observed colours of almost all galaxies. Mixtures with an in- 
termediate mass ratio between the old population and the 
burst produce also intermediate results in terms of indices 
and colours between the pure burst models (Fig. [5| and the 
CSP-S ones presented in Fig. [T] In fact, as shown in Fig. 
[S] a more moderate ratio between "old" and "burst" stars 
of 2:1 yields a decent match between the BC03 models and 
the observations. As for the Ma05 models, however, the fact 
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Figure 6. H/3 vs [MgFe]' diagnostic plot, showing that our galax- 
ies need to have metallicity not lower than ^ 40% solar in order 
to be consistent with the observed strength of metal sensitive ab- 
sorption lines. The four evolutionary tracks in colour are from 
BC03 models, colour coded and marked as in Fig. [5] 

that already SSPs are too red in optical-NIR colours results 
in no combination of "old" vs "burst" stars being able to 
reproduce the observed colours. 

We repeated the experiment with the 9:1 ratio, but 
changing the metallicity of the old component from Z — Zq 
fixed to a variable metallicity that matches the one of the 
"young" component. We tag these models as CSP-Y, where 
"Y" indicates that the metallicity of the old component 
is the same as in the young component. The results are 
shown in Fig. [9l BC03 models (Fig. [9] top) are able to cover 
most of the parameter space occupied by the data-points; 
however a significant number of galaxies are too blue in 
mo. 55 — mi. 40 with respect to models at given D4000n. More- 
over, the metallicities implied by colours are systematically 
lower (and to large extent inconsistent) with those implied 
by the absorption indices alone. As one can argue in the 
light of the results of pure SSP and 2:1 ratio models, less 
extreme fractions of the old component can provide a much 
better match to the data. As for the Ma05 models. Fig. [9] 
bottom shows that not even with this kind of mix it is possi- 
ble to match the distribution of mo. 55 — mi. 40 vs D4000n, the 
optical-NIR colours of the models being too red with respect 
to the data. We reiterate that any dust with the most com- 
monly adopted attenuation curves would only worsen this 
disagreement. 

In summary, adding a substantial fraction of old stel- 
lar populations to reduce the impact of TP-AGB stars pre- 
dicted by Ma05 in the NIR may indeed explain the lack of 
observed NIR features, while remaining consistent with the 
observed optical spectral indices. However, such a combi- 
nation fails to reproduce the observed optical-NIR colours 
at given II5a + H7A, D4000n and optical colours, no mat- 
ter if the metallicity of the "old" component is the same 



as the one of the "young" component or is fixed at the so- 
lar value (as it is probably more realistic for an old stellar 
population in a system of a few times 10^" M©). We cannot 
exclude that more complicated and ad hoc star formation 
histories and metallicity mixtures might be able to produce 
a closer match between observations and Ma05, but it ap- 
pears unlikely that such a fine tuning actually applies to all 
16 galaxies in the sample. On the other hand, BC03 mod- 
els appear to reproduce the observations much more easily 
by assuming SSPs or CSPs with even massive "old" compo- 
nents of variable sub-solar metallicity; in order for models 
with an "old" component at solar metallicity to match the 
observations, "old" fractions less than 0.9 are required in 
most cases: it is sufficient to lower such fraction to 0.6-0.7 
in order to obtain a generally good agreement. 

We further illustrate these points by directly comparing 
SSP and CSP model spectra with an actual spectrum of our 
PSB sample in Fig. [TO] where PSB J1006+1308 is used as 
a mere example which optimally serves to this goal, as it 
is broadly representative of the full sample. For both BC03 
(shown in the left-hand panel) and Ma05 (in the right-hand 
panel), we pick the best fitting models from each of the three 
grids of SSP, CSP-S and CSP-Y (assuming old:young=9:l). 
They are overplotted to the observed spectrum (solid black 
line), in blue for the SSP, in solid red for the CSP-S and 
dashed green for the CSP-Y. The "best fitting" models 
are those yielding the minimum o ver the optima l set o f 
optical absorption indices defined in iGallazzi et ahl (|2005l ) 
(D4000„, H5a + H7a, H/3, [MgFe]', [MgaFe]). Due to the 
coarse grid in metallicity, the use of fixed old fractions and 
the neglect of dust attenuation, we warn the reader not to 
over-interpret these fits in terms of derived physical param- 
eters. Nonetheless, the following key observations can be de- 
rived from Fig. 1101 i) All best fit models can reproduce the 
observed optical features extremely well (see insets). The 
slope of the optical continuum is also generally well repro- 
duced, although for BC03 the SSP fit is too blue in the opti- 
cal range, which might be due either to a small mismatch in 
metallicity because of the coarse grid or to small amounts of 
dust, which are not included in our simple models, ii) The 
NIR spectrum of the best fitting Ma05 SSP, with an age of 
0.6 Gyr and metallicity 2 times solar, has the typical NIR 
boosting due to the Ma05 implementation of the TP-AGB 
phase, clearly recognisable from the sharp features around 
1.41 and 1.77 nm. However, no evidence of such features 
appears in the observed spectrum; moreover, the observed 
spectrum is almost a factor of two dimmer than the Ma05 
SSP in the NIR. In contrast, the BC03 SSP provides a very 
good fit in the NIR both in terms of (lack of) features and of 
overall flux, in) The use of CSPs with a substantial fraction 
(90%) of 10-Gyr old stars yields a relative enhancement in 
the NIR over optical for BC03 models, for both CSP-S and 
CSP-Y. The plot suggests that a smaller old fraction and 
lower metallicity than the best fltting pure SSP might pro- 
vide a very close match to the NIR-to-optical ratio. CSPs 
for Ma05 (which by coincidence converge to the same so- 
lution for both CSP-S and CSP-Y, i.e. the solid red line 
and the dashed green one overlap) do actually yield much 
weaker NIR features and also lower NIR flux relative to op- 
tical. However, this stays significantly higher (« 50%) than 
observed. 
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Figure 8. Same as Figure[7l but now using CSP-S with a ratio of masses as formed in the fixed "old" vs variable "young" components 
of 2:1. This figure shows that a substantial underlying old stellar population is compatible with most of the observations if BC03 models 
are used, but produces inconsistent colours if the Ma05 are adopted instead. 



3.4 The anomalous galaxies PSB J1150 0107 and 
PSB J1206+0918 

PSB J1150-0107 and PSB J1206+0918 are characterised by 
extremely blue optical-NIR colours (mo. 55 — mi, 40 < 0.2) 
and intermediate/red optical colours (mo. 55 — mo. 70 0.3). 
When plotted on the colour-colour diagrams of Fig. [5] [T] E] 
and [51 they significantly depart from the sequences both of 
observed data-points and model tracks. None of them is able 
to reproduce such a combination of intermediate/red opti- 
cal colours and very blue optical-NIR ones: for the given 
mo. 55 —mi. 40, an offset of —0.1 to —0.2 mag in mo. 55 —mo. 70 
would be required to bring them back in agreement with the 
models, or, for the given mo.55 — mo.70, an offset of > 0.3 mag 
in mo.55— mi, 40. From the arrows plotted in the colour-colour 
diagrams, it is i mmediate to see that a standard effective 
attenuation a la ICharlot fc Fail l|200d ) is unable to move 
the models closer to these two objects: in fact, one would 
need an effective attenuation curve which is extremely flat, 
thus causing similar attenuation at 0.70 and 1.40 ^m. With 
standard dust grains, for which the absorption cross section 
decreases with wavelength, this can be only obtained by in- 



voking a peculiar geometry in which back scattering plays a 
major role. 

Significant residual star formation is witnessed in 
PSB J1150-0107 by EW(Ha)= -2.43 ± 0.24 and a 
marginal detection at 22 pm in the WISE all-sky survey. 
PSB J1206-I-0918, however, does not display any sign of star 
formation, with EW(Hq)= —0.65 ± 0.26 and no detection 
at 22 fim. This suggests that residual star formation or pro- 
cesses related to nuclear activity are not the cause of the 
anomaly common to these two galaxies. 

We note that the shape of the spectra of these two galax- 
ies is peculiar in two main aspects: i) the continuum appears 
to bend down, i.e. the slope becomes steeper, going from 
~ 0.55 to 0.8 ^im, whereas all other galaxies display an ap- 
proximate constant slope or even an up-bending shape in 
this region; it) a significant and featureless drop long- ward 
of « 0.5 fim is observed instead of the Fe and Mg absorption 
complexes, which is probably due to the low metallicity of 
these two galaxies, among the lowest of the entire sample. 
These two features appear to be responsible for the anoma- 
lous colours, by making the spectral slope between 0.7 and 
1.4 fim too steep relative to the one in the optical range 
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Figure 9. Same as Figure [7] for CSP-Y models (see text), composite SPs with the "old" component having the same metallicity as the 
"young" one, according to the corresponding track. 



0.55-0.70 /im. The fact that models fail to reproduce these 
colours may point to problems in the low metallicity regime, 
below « O.5Z0. 

We note that, according to the Ma05 models, the largest 
enhancement and the most dramatic features in the NIR 
are expected for Z < Zq. Instead these two low- metallicity 
galaxies exhibit bluer optical-NIR colours than even "TP- 
AGB light" models (e.g., BC03). 



3.5 Comparison with Kriek et al.(2010) 

As already mentioned in Sec. n iKriek et all (|2010l ') have 
obtained a detailed (rest-frame) UV-optical-NIR SED for 
a sample of 62 PSB galaxies observed in the NEWFIRM 
Medium Band Survey (NMBS) in the redshift range 0.7- 
2. Their composite SED is essentially a stack of the SED of 
their sample. They applied maximum likelihood fitting using 
both BC03 and Ma05 on different wavelength ranges (opti- 
cal only or optical plus NIR) and, similarly to the results we 
have presented in this work, they found that Ma05 cannot 
reproduce the blue optical-NIR colours observed, whereas 
BC03 provide an almost perfect match. 

In Fig. [11] we compare the SED published by 



iKriek et all l|201Ci ). represented by the red points with er- 
ror bars, with the stack of our 16 spectra. Both SEDs are 
normalised at 5500 A. The black line is the median of our 
16 spectra at each wavelength, while the shaded area dis- 
plays the 16-84 inter-percentile range. The similarity of the 
stacked SEDs from the two works is striking, typically con- 
sistent within 1-1.5 a. Our spectra are slightly bluer and 
steeper in the range 5000 A to 2 fim, which might de- 
rive from some relative bias in the respective sample se- 
lections., with our objects being somehow more extreme 
(e.g. in terms of burst fraction, see also Sec. In fact, 
iKriek et all HoiO) define their sample based on k-corrected 
rest-frame colour cuts which are hardly comparable to our 
spectroscopic criteria. Alternatively, this comparison might 
indicate that NIR-bright evolved (TP-AGB) stars have a 
slightly larger influence on the SED of PSB galaxies at 
z > 0.7 than at a » .2, co nsistently with expectations, e.g., 
of Melbourne et al.1 (|2012l ). Apart from this minor colour 
difference, our detailed spectroscopic study of individual 
PS B galaxies confirm s and reinforces the previous results 
bv lKriek et al] (|201G| ). based on stacking and medium band 
photometry and at a very different age of the Universe. 
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Figure 10. Example of how different models perform in reproducing the SED of a PSB galaxy, specifically PSB J1006+1308, whose 
observed spectrum is shown in black. In the left-hand panel models based on BC03 are overplotted, while models based on Ma05 are 
shown in the right-hand panel: the blue lines are used for pure SSP, red solid lines for CPS-S (90% 10-Gyr old component at solar 
metallicity) and green dashed lines for CSP-Y (90% 10-Gyr old component at the same metallicity as the young component). The models 
plotted here are chosen from each of the simple grids presented in the text and in Fig. [S] [7] and [9] as those that best fit the set of 
five optical stellar absorption features sensitive to age and metallicity defined in Gallazzi et al. (2005). All models can reproduce the 
spectral absorption features remarkably well (see insets) and agree on the presence of a young component approximately at the peak 
of the TP-AGB contribution. However their extrapolations to the NIR regime differ significantly. Most notably, while the inclusion of a 
substantial contribution from an old component can help reduce the strong NIR features predicted by Ma05 models, the NIR-to-optical 
ratio predicted by these models remains nevertheless too high. 



4 SUMMARY AND CONCLUDING REMARKS 

We have presented the first spectroscopic study of post- 
starburst galaxies covering both the optical (SDSS) and 
NIR (our new VLT-ISAAC observations) wavelength ranges, 
which was specifically designed to detect the two main sig- 
natures of the p resence of TP-AGB stars predicted by the 
iMarastonI l|2005l ) models in contrast to "TP-AGB light" 
models (such as BC03): i) the strong NIR absorption fea- 
tures caused by carbon composite molecular band-heads; 
and ii) the boosting of the NIR flux. To this goal, we have 
carefully selected a sample of 16 galaxies with post-burst 
ages corresponding to the peak contribution by TP-AGB 
stars to the NIR flux (0.5-1.5 Gyr), in a broad range of stel- 
lar metallicities (f» 0.2-2.5 Zq). These galaxies are at high 
enough redshift (0.15-0.25) to allow clean measurements 



of the NIR features at rest wavelengths 1.41 and 1.77 /im 
through the atmospheric H and K windows. Despite these 
ideal conditions for the detection of the characteristic signa- 
tures of TP-AGB stars predicted by the Ma05 models, the 
16 galaxies in our sample all display smooth, featureless NIR 
spectra and blue optical-NIR colours. Interestingly, these ob- 
served spectral properties are generally well reproduced us- 
ing the BC03 models assuming simple (i.e. single-age) stellar 
populations. 

We have investigated whether the observations could be 
reproduced with Ma05 models as well, by combining stellar 
populations of different ages. We find that diluting the light 
of a young (i.e. a few 100 Myr old) starburst with that of a 
massive old component allows the Ma05 models to reproduce 
the typical "E-f A" optical spectra of post-starburst galax- 
ies characterised by strong Balmer absorption lines and a 



© 2012 RAS, MNRAS OOO.mfSTI 



NIR spectroscopy of post-starburst galaxies 19 



.2 - 



0.6 - 



0.4 - 



0.2 - 




i Kriek et al. (20' 0)- 



0.0 r ii I 

0.125 0.250 



1.500 1.000 



2.000 '1.000 



Figure 11. Comparison between the stacked SED from our 
spectra and from the post-starburst sample from the NMBS of 
Kriek et al.(2010), shown as red points with error bars. The black 
line is the median spectrum of our 16 galaxies, whose spectra had 
been previously normalised at 5500 A. The shaded area shows the 
16—84 percentile range of the data. The two samples yield con- 
sistent results (although our sample appears marginally biased 
towards bluer SEDs relative to Kriek et al. 2010) and, most im- 
portantly, agree on the fact that there is no need for boosted NIR 
emission by TP-AGB as implied by Ma05. 



well developed D4000 break. At the same time, such mod- 
els exhibit very weak NIR features, in agreement with the 
observations, because TP-AGB stars from the young com- 
ponent have not developed yet, while the NIR emission from 
old stellar populations is smoother. However, these models 
still cannot account for the data, because the predicted ratio 
of NIR to optical flux is significantly larger than observed 
at the metallicity of our galaxies. Hence, the observed lack 
of deep NIR absorption features and the modest ratio of 
NIR to optical flux in our data cannot be simultaneously 
accounted for by a simple "light dilution" effect. We note 
that the observation of metal absorption lines at optical 
wavelength does not support the idea that the faint NIR 
continuum (and the consequent dilution and/or absence of 
strong NIR absorption features) could arise from extremely 
low metallicities. 

The Ma05 models predict too red optical-NIR colours 
for any choice of ratio between the old and the young compo- 
nent. On the contrary, models based on BC03 fail for "old" 
mass fractions as high as 90%, but produce colours that are 
consistent with (most of) the observations for "old" mass 
fractions up to 60-70%, which leaves enough space to accom- 
modate more reasonable SFHs than a (possibly unrealis tic) 
single burst. We note that the work of lPracv et al.l (|2012l ) in- 
dicates old fractions close to 90% for their seven PSB gala:x- 
ies, which would be inconsistent with our SEDs. However, 
this estimate might be biased by their assumption of fixed 
solar metallicity. Moreover, and probably more importantly, 
it is possible that our objects are somehow more extreme 
in terms of burst fraction, because, due to our redshift con- 
straints, we select objects close to the faintest limit of the 



spectroscopic sample of the SDSS. In this way we would se- 
lect the intrinsically most luminous PSB galaxies at 2 ~ 0.2, 
hence those with the largest burst fraction (i.e. with the low- 
est "old" fraction): in fact, for a given total mass, a higher 
burst fraction yields higher luminosity. The lack of strong 
colour gradients in most of our PSB galaxies lends support 
to this hypothesis, suggesting that the starburst substan- 
tially affected the f ull extent of t he ga laxy, rather than the 
centre only. Finally, iPracv et al.l (|2012f ) point out that their 
sample is less extreme than typical samples of PSB galaxies 
in the literature due to the small physical size of the spectro- 
scopic fiber (few hundred parsecs), which bias their selection 
towards central rather than extended starbursts. In conclu- 
sion, the consistency of our SEDs with BC03-based CSPs 
and "old" mass fractions up to 60-70% does not raise any 
problem for the physical interpretation of our PSB galaxies 
using the BC03 models. 

We also find that attenuation by interstellar dust, using 
standard prescriptions, can only worsen the disagreement 
(when present) between spectral evolution models and the 
observed optical-NIR colours of post-starburst galaxies in 
our sample. This is particularly true for the Ma05 models, 
which predict systematically redder colours than observed. 
We note that heavier-than-implemented attenuation by cir- 
cumstellar dust of the most luminous C stars about to eject 
their envelopes could potentially reduce the contribution by 
these stars to the integrated spectrum and perhaps improve 
the agreement of the Ma05 models with observations!^ A 
more accurate theoretical and empirical assessment of the 
effects of circumstellar dust absorption is definitely key to 
make progress in the SED modellin g of the stellar popula- 



tion affected by TP-AGB stars (e.g . Srinivasan et al.l |2009| : 



iMelbourne et al.ll2012i : iMeidt et ani2012l ') 



The results presented i n this work confirm the previous 
claim bv lKriek etHI l|2010l ) that BC03, and more generally 
"TP-AGB light" models, provide a better agreement with 
the spectral properties of galaxies hosting a stellar popula- 
tion at the peak of the TP-AGB phase than Ma05 models. 
The important novelty of our work is to base these conclu- 
sions on a higher-resolution, high-quality spectrophotomet- 
ric dataset and an optimal spectroscopic sample selection. 

Somehow in contrast with these studies and our work 
in particular, iLvubenova et all (|^012) have detected the 
1.77 fj,m feature in a couple of Magellanic Cloud's globu- 
lar clusters and found them in agreement with Ma05 pre- 
dictions. This result calls for more observations in order to 
obtain a representative sample of globular clusters over a 
broader range of metallicities (the two in which the detec- 
tion was obtained are both at metallicity between 0.2 and 
0.6 solar) and beat down stochasticity effects: in fact, a sin- 
gle (TP-)AGB star can produce a substantial fraction of the 
total cluster NIR luminosity in the mass range explored. 



The Ma05 mod els rely on period-averaged s pectra of TP-AGB 
stars assembled bv lLancon &: Mouhcinj 1I2OO2I ). who recommend 
reddening these spectra before use in population synthesis mod- 
els (see their section 3.3). Whether this or a different attenuation 
has been included cannot be asserted from the description of the 
TP-AGB spectra in section 3.3.2 of Ma05. In the BC03 models, 
the circumstellar attenuation of TP-AGB stars experiencing su- 
perwinds is accounted for in an empirical way (see section 2.2.4 
of BC03 for detail). 
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One further possi bility to reduce the apparent disagreement 
between the L vubenova et all l|2012h results and ours might 
reside in a technical aspect of their challenging observations: 
the addition of targeted AGB stars outside of the central re- 
gion to integrate the central spectrum and obtain full spatial 
coverage of the globular clusters might lead to AGB stars 
being over-represented with respect to the stars that con- 
tribute the bulk of the diffuse emission. 

In the light of this complex and partly controver- 
sial observational framework, the apparent agreement of 
the BC03 models (of either simple or composite stellar 
populations) with the observations presented in this pa- 
per should not be necessarily interpreted as an indica- 
tion that the description of the TP-AGB phase in these 
models is adequate. New models incorporating more re- 
cent developme nts in the theory and observation of TP- 
AGB stars (e.g. iMarigo et al.l l2008l: ISrinivasan et ahllioogl : 



iGirardi et a!]|2010l : lMelbourne et al.H2012D are being assem- 
bled (Chariot & Bruzual, in preparation). By providing a 
better account of observations, new models of this type will 
allow deeper insight into the physical properties of galax- 
ies from spectral interpretations, especially at high redshift, 
where a large number of galaxies is going through the evolu- 
tionary phase affected by TP-AGB stars. In this perspective, 
the present work provides key and very stringent constraints 
for the new generation of stellar population synthesis mod- 
els. 
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